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Abstract 

This paper gives examples of explicit arbitrage-free term structure 
models with Levy jumps via state price density approach. By gen- 
eralizing quadratic Gaussian models, it is found that the probability 
density function of a Levy process is a "natural" scale for the process 
to be the state variable of a market. 
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1 Introduction 

1.1 Literature review 

In the classical Black-Scholes economy, the exponential of a Brownian motion 
(with drift), i.e., 

S t = S exp(aB t + bt), 

is used to model stock prices. When one takes jumps into account, extending 
the Black-Scholes scale naturally leads to the modelling economic factors by 
exponential of a Levy process, say: 

X t = X eyipZ t , 
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for Scientific Research (No. 18540146) from the Japan Society for Promotion of Sciences. 
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where Z is a Levy process. Such a model is often called "exponential-Levy" 
type, and has been widely used in financial modelling since Merton's jump 
diffusion model [Mcr76j appeared. The variance gamma model by Madan and 
his co-authors (see e.g. [MCC98j ). and the hyperbolic model by Eberlein and 
his authors (see e.g. |EP02j ). are two major representatives, but there are 
still many others including [BN98 ] and [BNLOlj . For further references, see 
|Miy01| or [CTU4l Chapter 8-11]. 

In the context of interest rate models, the exponential scale is also com- 
mon (though the scale is in general not time- homogeneous any more); e.g., 
bond prices of Vasicek's model |Vas77j or Gaussian HJM [H.TM92] , LIBOR 
rates of the BGM's LIBOR market model |BGM97| . etc. There are also 
exponential-Levy type interest rate models with jumps; e.g. Shirakawa's ex- 
tension of Gaussian HJM [Shi91j, Chiarella-Nikitopoulos [CNS03] . Eberlein- 
Raible's Levy term structure models |ER99] . or Albeverio-Lytvynov-Mahnig 
|ALM04j . etc, etc (see also section [2.2.31 in the present paper). 

There are two major substitutes in the continuous-path cases; the affine 
class and the quadratic class. Roughly speaking, we say the model is "(exponential- 
quadratic" (w.r.t. a state process Z) if the price at time t of the bond with 
maturity T is given by 

= exp{(A 2 (t,T)Z t ,Z t ) + (A 1 (t,T),Z t ) + A (t,T)}, 

where A%{t, T) G R n <8> R n is a non-zero symmetric matrix, A\(t, T) is in R n , 
and Ao(t,T) G R. If Ai = then the model is called "(exponential-)affine". 

Affine class w.r.t. a Levy process is nothing but the above-mentioned 
exponential-Levy class. When we refer to affine class, the state process is 
usually a Markov process whose state space is a convex cone in R n (and hence 
symmetric processes are excluded). The class contains Cox- Ingersoll- Ross 
model [Cl R85j and Duffie-Kan's multi-dimensional generalization [DK96J. 
The jump-diffusion extension of the affine class was fully done by Duffie- 
Filipovic-Schachermayer |DFS03j by generalizing Kawazu-Watanabe theo- 
rem [KW71] . It is also shown by Filipovic-Teichmann [FT04] that the " finite- 
dimensional realization" is limited to the affine class. 

The quadratic class, which can be embedded into the (degenerate) affine 
class, is gaining much popularity among economists and practitioners mainly 
due to its tractability, (see e.g. AI)G()2.) or partially explained from a 
mathematical background; it has a rich mathematical structure (sec [AH06J). 
Contrary to the affine case, however, no jumps are allowed in the quadratic 
class as Chen-Filipovic-Poor showed in [CFP04J within Markovian frame- 
work (though Levendorskh |Lev05| gives a jump-type extension of quadratic 
models "in the direction of affine class"). 
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The main aim of the present paper is to give a description of a Levy 
extension of the quadratic models in a totally different way from |CFP04j 
and |Lev05] . 



1.2 Our result 

The above story of the quadratic term structure models tells us that exponential- 
quadratic functions are "natural" scales (at least in modelling interest rates) 
for Brownian motions but not for pure-jump Levy processes. Then what 
scale is natural for a Levy process? 

To answer this question we start from the observation that within the 
exponential-quadratic scale lies the probability density of the Gaussian dis- 
tribution. This observation naturally leads to asking if its probability density 
is a natural scale for a general Levy process or not. To be more precise, the 
question is: given a Levy process Z, the market consists of 

Pj = p(A(t, T), Z t ) (or something like this) 

where p(t,x) = P(Z t G dx\Z = 0)/dx can be consistent? or in particular, 
arbitrage-free? 

An answer to this question is given as Theorem 13.11 It says that if the 
instantaneous forward rate f(t,T) := —dx^ogP^ is given by 



-d T logp(X T + T -t, Z, 



for some continuous A : [0, oo) — > [0,oo), then the market is arbitrage-free 
(Theorem 13.11) . 

To construct arbitrage-free interest rate model, we rely on so-called state 
price density, or pricing kernel approach [BMOll pp371-], which was initi- 
ated by Constantinides [Con92j and later developed by Rogers |Rog97| and 



Hughston and his co-authors [FH9^ |BH04[IHR05] (see also the textbook by 
Hunt and Kennedy |HK04 ]). 



1.3 Organization of the present paper 

We will start from a brief survey of the approach (section l2.ip . and then give 
two important classes; Gaussian (section 12. 51) and quadratic Gaussian (sec- 
tion ^. 2. 21 . A direct jump-type extension which we call generalized Shirakawa 
model will be given in section 12.2.31 In section [3TT1 Theorem 13.11 and its proof 
will be given. Examples based on our new framework will be presented in 
section 13.31 
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2 The state price density approach to inter- 
est rate modelling 



2.1 Review 

In principle, a strictly positive process {nt} is a state price density with 
respect to a market on a filtered probability space (O, J 7 , P, {^t}) if for any 
asset indexed by i G / that generates {D l s } cash flows in the future, its 
market price at time t is given by 

POO 

Sl = n^E p [ K s dDi\F t ], (2.1) 
Jt+ 

or for any T > t, 

Sl = 7ir 1 E p [n T S l T + f K s dDi\F t \. 
Jt+ 

In other words, HT/^t ( x probability density with respect to P) gives (ran- 
dom) discount factor of a (random) cash flow at time T. In particular, if we 
denote by Pj the market value at time t of zero-coupon bond with maturity 
T, we have 

P? = ^ 1 E F [n T \F t }. (2.2) 

In the famous text book by Duffie |Duf01j . an arbitrage- free market is 
characterized by the existence of a state price density. This duality, which 
comes from the Hahn-Banach theorem and its variants, is proven for fully dis- 
crete markets (finite dimensional cases) and for Brownian markets (Brownian 
filtration cases). In more general cases one can extend the duality though 
she/he needs to be careful about the exact meaning of the arbitrage. It de- 
pends on what kind of assets are traded and what kind of trading strategies 
are admissible in the market. To determine how far we can extend is, how- 
ever, out of the scope of the present papeiQ. Here we just assume that the 
space of the value process of trading strategies is orthogonal to it in the L 2 
space of stochastic processes, meaning that we presume that the existence of 
a state price density implies the market is arbitrage-free (though we leave the 
problem of market completeness behind). This assumption is robust because 
at least it is fulfilled if we admit simple strategieJ^ only. 

1 Extensive studies in this direction are found in DS94, DS99, DS06 , and for the cases 
of bond markets [BMKR97] . |BKR97j . and [DD04] to name a few. 

2 by which we mean those strategies which remain constant for a short time interval. 
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From a perspective of modeling term structure of interest rates, the for- 
mula (12. 2p says that, given a filtration, each strictly positive process n gen- 
erates an arbitrage- free interest rate model. On the basis of this observation, 
we will construct arbitrage-free interest rate model. 



2.2 Examples 

In practice, an explicit formula for the bond price (and hence interest rates) 
is desirable; it gives a parameterization of the entire term structure, and 
it becomes applicable to, say, duration-based hedging. The more we have 
explicit formulas for present values S l for the cash flow D l through ( 12. ip the 
better the model becomes. 

Below we give two classical Term Structure Models (TSMs for short) 
which exhibit no jumps. 



2.2.1 Gaussian TSMs 

Let W ( be a standard Brownian motion taking values in H d and h(t, s) be 
an element of H <g> H where 

H = {h : [0, oo) — > R d , absolutely continuous s.t. h <E Lf oc [Q, oo)}. 

Given an initial data T > Pq, define 

7r t = P *exp |jT (h s (t,s), dW s ) Rd - ~J \h s (t,s)\ 2 ds 
Then, by an easy manipulation we have 



(2.3) 



i r* 



(2.4) 

{\h s (T,s)\ 2 -\h s (t,s)f)ds 



where h s stands for the partial derivative with respect to the latter variable. 

Note that this class covers so-called Gaussian Heath- Jarrow-Morton mod- 
els (see e.g. |MR97l 319-]). In fact, we have 



f(t,T) :=-d T \ogP? 
= /(0, T) + J (-h TjS (T, s), dW s ) Rd + J (h T;S (T, s), h s (T, s)) Rd ds ^ 5 ^ 

= /(0,T)+ / (-h T , s {T,s),dW s -h s (T,s)ds) nd . 
Jo 
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The last expression shows that {f(t,T)} t <T is a martingale under the so- 
called forward measure P T defined by 

P T 

alP T /dP- 1 



E[P t T ] 

= exp|-^ (h s (T,s),dW s ) nd -~^ \h s (T,s)\ 2 ds\. 

2.2.2 Quadratic Gaussian TSMs 

Let A be a continuous map on R + taking values in the set of all positive 
definite d x ci-symmetric matrices, and k : R + — > R be a continuous map. 
Define 

Tc t (x) = exp {-(A t x,x) R d + k t } , (x G R d ). 
Then we have, for a d- dimensional Wiener process W starting from the origin, 

Pj = E[7i T (W T )\a(W s] s < *)]/7r t (W t ) 
= {det(2(T-t)A T + /)}- 1 / 2 - 

• exp (-({A T -At- 2(T - t)(2(T - t)A T + J)- 1 } W t , W t ) + (k T - k t )) , 

(2.6) 

where / is the unit matrix. 

For the derivation of (12.61) . see Appendix. 

Remark 2.1. There are several ways to introduce QTSMs. The above is 
just an illustration. For details see e.g. [CFP04j and [A"H06]. 

2.2.3 Generalized Shirakawa TSMs 

Let p be a stationary Poisson point process on a measurable space (E,Be)- 
Then its counting measure defined by 

Af p ((s : t],A) := G (8,t];p(u) G A}; (s < t, A G B E ) 

is a stationary Poisson random measure; i.e. for mutually disjoint Bi, B n G 
Bh +x e, the random variables J\f p (BjYs are mutually independent and Poisson 
distributed, and for A G Be and s < t, 

e -{t-8)v{A) 

P{Af p {(s,t],A) =n) = {(t - s)v{A)Y : , 

n! 

where u is a a- finite measure on (E, Be)- For a detailed instruction, see e.g. 
jIW89] . 
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Let 5 be a Borel function from R + x R + x E to R such that (the equiv- 
alence class of) S(t, -, •) is in n p >iL p (dt ® z/) and sup p>1 \\S\\ P LP < oo. 
For this 5 and /i G H <8> H , define 

Z t = h s (t,s)dW s — - I \h s (t, s)\ 2 ds + / / 8(t,s,x)J\f p (ds,dx), 
Jo 2 Jo Jg J^ 

where if, /i s , and W are as in Example 12.2.11 and W is independent of p (or 
equivalent ly, Af p ). 

By the assumption on h and 5, we have 

£[e Zf ] = exp jjf y (e^-.*) - l)i/(da;)dsj 

= exp < / / {5(i, s, a:)} p j/(da;)ds > < oo. 
[ p=1 Jo Je J 

Given an initial data T i— > Pg 1 ", define 

vr, := P*e^/-B[e^]. 
Then, denoting the bond price in Gaussian TSM (12. 4p by p t T > Gauss ; we have 
= E[n T \rt W ]/n t 
= p t T ' GdUfah eX p | J J (J(T, s,x) — 6(t, s, x)) J\f p (ds, dx)^ 



• exp 



o JE 



Here the filtration {J-f ,W } is the one generated by p and W. 

If we further assume that 6(t, •, •) is differentiable in t and if its derivative 
d t 5{t, -, •) is, say, uniformly bounded and if v is finite, then denoting the 
forward rate in Gaussian TSM (12.51) by f Gauss (t,T), we have 

f(t,T):=-d T \ogPj 

/■* r (2 7) 

= / Gauss (t, T) + / / *r(T, s, x) {A/" p (cfe, da;) - e S{T ' s ' x) v(dx)ds} . 
Jo Je 

This expression can be regarded as an extension of Shirakawa's model iSliif) 1 i , 
and as a special case of jump-diffusion or fairly general semi-martingale mod- 
els, e.g. by |BKR97j . 
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Remark 2.2. If E is a finite set, then for each t > 0, uj G Vt we define the 
Poisson integral of 5 as a random finite sum by 

I 5(x)N(t, dx)(u) = S(x)N(t, {x})(cu), 

and in this case (12.71) becomes Shirakawa's model. 

The generalized Shirakawa model we have presented is basically a exponential- 
Levy model, and therefore we cannot use those Levy process without finite 
moments. For example, we miss the symmetric a-stable processes, which is 
characterized by (constant times) |£| Q (a G (0,2)) as its Levy symbol, i.e.; 
the Levy process Z with the property 

E[exp{i(£, Z t - Z s )\ = exp{-(t - s)6\£\ a }, 9 > 0, £ G R d , < s < t. (2.8) 

In fact its a-th moment explode. To construct an interest rate model driven 
by stable processes, we might use an approach given in the next section. 



3 Levy Density TSMs: A Generalization of 
QTSMs 

3.1 Main result 

Let Z be a Levy process in R d starting from the origin, adapted to a given 
filtration {J-'t}- We assume that it has the probability density p(t,x) with 
respect to the Lebesgue measure of R d : 

P(Z t G dx) = p(t,x)dx. (3.1) 

Here we assume that 

P(p{s,Z t ) > 0, Vs > 0,Vt > 0) = 1, (3.2) 

and 

p(t, -)p(s, •) is integrable w.r.t. the Lebesgue measure for all t,s > 0. 

(3.3) 

Then we have the following. 

Theorem 3.1. Under A3, jjj) and \3. 'J\) , the term structure model given by 

P?=p(\ T + T-t,Z t )/p{\t,Z t ), 0<t<T<oo, (3.4) 

where A : [0, oo) — > [0, oo) is a continuous function, is arbitrage- free in the 
sense of ( Iff.ffl) by putting n t = p(\t, Z t ). 
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Proof. By the Markov property of Z, we have 

E\p{\ T ,Z T )\T t ] = E\p(\ T ,Z T )\a(Z t )] 

= J ^p(\ T ,x + z t )P(z T _ t e dx) 

p(Xt, x + Z t )p(T — t, x) dx. 



R d 



Denote Levy symbol of Z by ip: 

JTp(At, •)](£)= / e^p(X T ,x)dx = e' x ^\ (^R d ). 

jR d 



Then we have 

p(T -t,x)= ^*[e-( T -^»](x) (= (2n)- d J e -HM e -{T-me) d A . 

Thus we can rewrite ( 13.51) as 

E\p(\ T ,Z T )\F t ] = (p(A r ,- + ^),^[e-( T -^)]) L2 . 
Since T o JF* = id on the space of density functions, we have 

<p(A r , • + Ztl^ie-V'^Dv = (2n)- d (Fip(\ T , ■ + Z t )] l e-P , - t W>) L .. 
Observing that 

JTp(Ar, ■ + Z t )](x) = e-^ z ^[p(X T , -)](x), (x G R d ), 

we have 

S[p(A T ,Z T )|^] = (27r)- d (e- J <-^e- A ^,e-( T -^-)) L2 

= (2rr)- d [ e^ x ^e-^ +T ~ l ^dx 

Jn* (3.6) 
= >F *[ e -(A T +r-t)W]( Zt ) 

= p(A T + T 

Hence we have 

£[p(A r , Z r )|^ f ]/p(A t , Z t ) = right-hand-side of (Ell). 
This proves the assertion. □ 
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3.2 Minor extensions 



3.2.1 Translation 

The assumption that Z is starting from is just a convention. The result is 
stable under the change of starting point. More precisely, for any z G R d , 
substituting Z t + z for Z t in (13.21) and (13.41) does not cause any problem. 

3.2.2 Product 

Let Z X ,...,Z be mutually independent Levy processes starting from the 
origin, A' : R + — * R + be continuous maps, and pi(t, x), I — 1, . . . , k be their 
probability density functions satisfying ( 13. 2\\ and ( 13. 3D . Then from Theorem 
13.11 it is clear that the bond market given by 

k 

Pi = \{pi^t + T -t,Z l t )/ Pl (t,Z l t ) 
i=i 

is arbitrage-free. In fact, one can take Y\i=iPi(^j ^t) to be a state price 
density. 



3.3 Examples of LDTSMs 

3.3.1 Quadratic Gaussian TSMs as LDTSMs 

Let Z be a d-dimensional Gaussian process such that Z t+ At — Z t ~ N(0, £At), 
where E is a positive definite matrix. Then its probability density is given 
by 

p(t,x) := (27rt)^ /2 (detS)- 1/2 exp{-^(S- 1 x,x)}. 

Let Xt : R+ — > R+ be an increasing continuous function. Using the formula 
( 13.41) . we have an LDTSM generated by Z as 



Pj = p(X T + T-t,Z t )/p(X t ,Z- l 



A, \ d/2 1 (3-7) 



Xr + T _ tJ ex P{-2^ Ar + T " £ ) " A * K E 
Since W t := S _1 / 2 Z is a Brownian motion, (13.71) is also represented as 

/A \ d/2 1 
Pi = ( Ar + r_ J exp{--{(A T + T - t)- 1 - A^jlW,! 2 }. (3.8) 

Note that := |W 4 | 2 is a <i-dimensional squared Bessel process, and so (13.8j) 
is a 1-dimensional affine TSM. 
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Set 

A t :=— I and k t := -\ log(27rA t ) d . 
zAj z 

Then we have 

p(X t , Z t ) = exp {-(A t W t , W,) Rd + kt} . 

Therefore, we could also obtain (13.81) using the formula (12.61) . 

It should be noted that a more general QTSM like (12.61) is within our 
LDTSMs; it is shown via the "product" construction of section 13. 2. 2[ 



3.3.2 Cauchy TSMs 

Let Z be a Cauchy process in R d starting from 0, whose probability density 
is given by 

/ v / , x r((d+i)/2)ot , s 

P(i , x) := P(Z« 6 = (r{(w)2 l ; |l _ ) ; | ) 2}) , J+1 , /2 . (3.9) 

where # > and 7 G R d . Note that Cauchy processes are strictly stable (or 
self-similar) [f| process with parameter 1; namely, 

zT ct = czT t , c>0,t>0. 
This property is seen from their Levy symbol: actually we have 

Apparently p(t, x) in (13.91) satisfies the assumptions (13.21) and (13.31) . Thus 
we can obtain a closed form expression for an arbitrage-free TSM driven by 
a Cauchy process as follows 



pT _\ T + T-tf 9 2 \ 2 t + \Z t -t^ 



t 



x t \e\\ T + T-ty + \z t -t^ 



(d+l)/2 



3.3.3 Gamma TSMs 

Let Z be a one dimensional gamma process with parameters a, b > 0, so that 
each Z t has density 

hat 



3 About the stable distributions, |Sat99j is a good reference. 
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for fixed x > 0. It is easy to check that (13. 2p and (13.31) are satisfied. The 
market value at time t is given by 

pT = JjXt —At +T—t r(aA0 a(AT - W) 

1 (a(\ T + 1 —t)) 

Here we may think Zq > (see section [3.2. lj) . 

The expression Pf jumps only upwards, which is unrealistic for the bond 
price movement. But as we remarked in section 13.2.21 we can consider Pj in 
dsn factor of the bond price. 

Remark 3.1. We can further generalize the Gamma TSMs to include CME 
(convolutions of mixtures of exponential distributions) and others appearing 
in |Yam06] . 



Appendix: A derivation of ( 12.61) 



R d 



We let kt = for the moment. By the Markov property of W, we have 
E[n T (W T )\a(W s ;s < t)] = E[ir T (W T )\a(W t )] 

tt t (.t)P(W t G dx\W t ) 

= [ ir T (x + W t )P(W T _ t edx) 

jRd (3.11) 

= (m^)) L Mx + Wt) exp ("2(r=l)) dx 

The last equality comes from the following: 

(A T (x + W t ), (x + W t )) Rd = (A r W t , W t > + 2(A T x, W t ) + (A T x, x) R *. 

Let Ut be an orthogonal matrix such that At = U T diag[e T , e T ]Ur- 
Here, of course, e T , are eigenvalues of A T . Then we have 

2(A T x,W t ) + {A T x } x)^ d 
= 2(diag[e^, e T ]U T x, Ef r W t ) R d + (diag[e^, e T ]U T x, U T x)^ d . 

Thus, putting y = {y l t , ...,yf) = U t x, 

I X I ^ 

2(A T a;,W t ) + (A T x,x) Rd + 



2(T - t) 

2(diag[4, 4]l/, ^W() R ^ + (diag[4, e T ]y, y) Rd + 



\y 



2 



2{T-ty 
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and denoting q = Cj(t, T) := e l T + 2 <T-t) ' 
= ^(Q(yl) 2 + 2^(^W t ) t ) 

2=1 

t=l L ^ C i ) i=1 



We claim here that 



It = 2(T _ t) ((2(T _ t) ^ + j)-i Wt) Wi); (3 12) 



1=1 



where / is the unit matrix. 

Assuming (I3.12p for the moment, we have 

(the last expression of (13.111) ) 

± \ -({A T -2(T-t)(2(T-t)A T +/)- 1 }w t ,W t )_ 



2vr(T - 1) 

1 \ d/2 

1 » e _({A T -2(T-t)(2(T-i)A T +/)- 1 }w t ,W t > / _ _ Q )-l/2 



2(T - 1) 

= {det(2(T - t)A T + j)}-V2 e -<{^-2(T- t )(2(T- t )^ + /)-}w tl w t >_ 

Thus we have (12.61) . 

Now we prove (13.121) . First observation is the following: 

diag[ci, c d ] = diagfe^, 4] + 2 (T-t) 1 

= U T A T m + — - — -/ = Ut(At + -t^- — r/)E#. 

T 2{T-t) v 2{T-t) ' T 

We then notice that 

diag[ I,...,I ] = f/T( ^ + _L_ /) -^. 
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Therefore, 



-A (U T W t ) 2 1 1 

£ 1 T tJf = (diag[-, -}U T W t , U T W t ) 

= (t^diag[-,...,-]C/ T W t ,W t ) 

= ((^T + ^ 7 y/)- 1 W,,W,) 

= 2(T - t)((2(T - t)A T + iy'Wt, W t ). 



This completes the proof. 
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